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In Hofava-Lifshitz gravity, regularity of a solution requires smoothness of not only the spacetime 
geometry but also the foliation. As a result, the regularity condition at the center of a star is more 
restrictive than in general relativity. Assuming that the energy density is a piecewise-continuous, 
non-negative function of the pressure and that the pressure at the center is positive, we prove that 
the momentum conservation law is incompatible with the regularity at the center for any spherically- 
symmetric, static configurations. The proof is totally insensitive to the structure of higher spatial 
curvature terms and, thus, holds for any values of the dynamical critical exponent z. Therefore, we 
conclude that a spherically-symmetric star should include a time-dependent region near the center. 
This supports the picture that Hofava-Lifshitz gravity does not recover general relativity at low 
energy but can instead mimic general relativity plus cold dark matter: the "cold dark matter" 
accretes toward a star and thus makes the stellar center dynamical. We also comment on the 
condition under which linear instability of the scalar graviton does not show up. 



I. INTRODUCTION 

Recently, a power-counting renormalizable gravity the- 
ory was proposed by Hof ava [l|, • The essential reason 
for the power-counting rcnormalizability is that in the 
ultraviolet (UV), the theory exhibits the Lifshitz-type 
anisotropic scaling 

t b z t, x bx, (1) 

with the dynamical critical exponent z > 3. Because of 
the Lifshitz scaling, this theory is often called Hofava- 
Lifshitz gravity. Although power-counting renormaliz- 
ability does not necessarily imply renormalizability, there 
is a good possibility that the theory is unitary and renor- 
malizable. Evidently, this is one of the driving forces be- 
hind recent enthusiasms for research on various aspects 
of the theory such as cosmology d, Q , black holes @ , 
the solar system dynamics @, and so on 

In Hofava-Lifshitz gravity a black hole is not "black" 
for high-energy probes even classically while it is indeed 
"black" for low-energy probes. The Lifshitz scaling im- 
plies that the dispersion relation is of the form uj 2 cx k 2z 
and the group velocity v g cx fc z_1 is superluminal in the 
UV. This means that, although there are "black hole" 
solutions in this theory, waves with sufficiently high fre- 
quencies can probe the region deep inside the gravita- 
tional radius and the central singularity can be seen. On 
the other hand, waves with low frequencies follow the or- 
dinary dispersion relation u 2 cx k 2 and thus cannot probe 
the region inside the gravitational radius. In this sense 
a black hole horizon is only an emergent concept in the 
infrared (IR). 

In order to investigate what happens inside a black 
hole, especially near the singularity and to see if the clas- 
sical singularity can be resolved by quantum effects, we 
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need to formulate a regular initial condition and to evolve 
it dynamically towards gravitational collapse. 

The purpose of this paper is, as a first step towards this 
outstanding problem, to initiate investigation of stellar 
solutions in Hofava-Lifshitz gravity. This is important 
also for constraining the theory from astrophysical ob- 
servations. A natural strategy would be to try to con- 
struct a globally static solution regular at the center. 
Surprisingly, we find a no-go result: we prove that a 
spherically-symmetric, globally static solution can not be 
constructed in Hofava-Lifshitz gravity under the assump- 
tion that the energy density is a piecewise-continuous, 
non-negative function of the pressure and that the pres- 
sure at the center is positive. Therefore, under the as- 
sumption in the matter sector, a non-singular stellar solu- 
tion should include a time-dependent region, presumably 
near the center. This supports the picture that Hofava- 
Lifshitz gravity does not recover general relativity at low 
energy but can instead mimic general relativity plus cold 
dark matter: the "cold dark matter" accretes toward a 
star and thus makes the stellar center dynamical. 

The rest of this paper is organized as follows. In sec- 
tion [H] we shall describe the definition, properties and 
problems of Hofava-Lifshitz theory. In section Mil we 
present the main result of this paper, i.e. the proof 
of non-existence of spherically-symmetric, globally-static 
solution regular at the center under a set of reasonable 
assumptions in the matter sector. Section [IV] is devoted 
to a summary of this paper and discussion. In appendix 
lAl we present asymptotically-flat vacuum solutions with 
A = 1. Those solutions presented in appendix [X] are 
used in appendix [B] where we show yet another no-go 
result for globally-static star solutions with A = 1. In 
appendix [C] we derive the condition (|24[) under which 
linear instability of the scalar graviton does not show up. 
In appendix [D] we consider a simple stellar solution in 
general relativity in Painleve-Gullstrand-like coordinate 
system to see how significant the regularity of time folia- 
tion is and to foresee what happens near a stellar center 
in Hofava-Lifshitz gravity. 
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II. PROPERTIES OF HORAVA-LIFSHITZ 
GRAVITY 

In this section we shall describe the definition, proper- 
ties and problems of Hof ava-Lifshitz gravity. 



A. Definition and basic feature of Hofava-Lifshitz 
gravity 

Hofava-Lifshitz theory is fundamentally non- 
relativistic since the Lifshitz scaling treats time 
and space differently. To be more precise, the funda- 
mental symmetry of the theory is invariance under the 
foliation-preserving diffeomorphism: 



t^t(t), x i -> x\t,x). 



(2) 



The foliation-preserving diffeomorphism consists of the 
space-independent time reparametrization and the time- 
dependent spatial diffeomorphism. This means that the 
foliation of the spacetime manifold by constant-time hy- 
persurfaces is not a gauge but has a physical meaning. 

Basic quantities in the gravity sector are the lapse 
function N(t), the shift vector N l (t,x) and the three- 
dimensional spatial metric gij(t,x). These variables can 
be combined to form a four-dimensional metric in the 
Arnowitt-Deser-Misner (ADM) form: 

ds 2 = -N 2 dt 2 + g %3 (dx l + N l dt){dx ] + N j dt). (3) 

Since the lapse function is roughly speaking a gauge 
freedom associated with the space-independent time 
reparametrization, it is rather natural to restrict the 
lapse function to be space- independent. This condition, 
called the projectability condition, is not only natural 
but also mandatory, as pointed out in Hofava's original 
paper Indeed, if we abandoned the projectability 
condition then we would face phenomenological obsta- 
cles [§[ and theoretical inconsistencies @. On the other 
hand, with the projectability condition (and without the 
detailed balance condition), the theory is free from those 
problems [3] . (See the discussion about strong coupling in 
the next subsection.) Therefore, throughout this paper, 
we impose the projectability condition and demand that 
the lapse function be space-independent. The Hamilto- 
nian constraint is, as a result, not a local equation satis- 
fied at each spatial point but an equation integrated over 
a whole space. 

Under the infinitesimal transformation 

St = f(t), 5x l = C(t,x), (4) 

gij , N* and N transform as 

$9ij = fdtgij + C-c9ij 

5N* = d t {N l f) + d t C +C C N\ 

5(Ni) = dt(Nif) + 9ij dtC j + C ( N h 

SN = d t (Nf). (5) 



Thus, N remains independent of spatial coordinates after 
the transformation. In the IR, where dt and dx % have the 
same scaling dimension, it makes sense to assemble g^, 
N 1 and N into a 4-dimensional metric in the ADM form 
([3]). The action is 



(6) 



where 



K, 



2N 



J dtdx z N^g{K l] K t] - XK 2 

+A + R + L Z>1 ), (7) 



(d t9ij - DiNj - D s N t ), l< .'/'*' A',,- (8) 



Di is the covariant derivative compatible with g^ , A is a 
cosmological constant, R is the Ricci scalar of g^, L z> \ 
represents higher spatial curvature terms and I m is the 
matter action. Here, we have rescaled the time coordi- 
nate so that the coefficients of K l] Kij and R agree. Note 
that not only the gravitational action I g but also the 
matter action I m should be invariant under the foliation- 
preserving diffeomorphism. 

By variation of the action with respect to N(t), we 
obtain the Hamiltonian constraint 



H„± + H m ± — 0, 



(9) 



where 



J dx 3 y/gn g ±, 



5Ig 

6N 

■W= I ''■•■ i vT/T I . no) 



and 



M 



^(Kv pij -A-R-L z>1 ), 



Here, pij and are defined as 



Pij =Kij- XKg 



(11) 
(12) 



and 



n^dx» = -Ndt, = ^(dt~ ( I •'! ) 



Variation with respect to N z (t, x) leads to the momentum 
constraint 



where 



Ho 



Hmi 



y/gSN* 
1 5I m 



(14) 



(15) 
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Note that the momentum constraint is determined solely 
by the kinetic terms and thus is totally insensitive to the 
structure of higher spatial curvature terms. In particular, 
for A = 1 the momentum constraint agrees with that in 
general relativity. 

As in general relativity, the gravitational action can be 
written as the sum of kinetic terms and constraints up to 
boundary terms: 



dtdx 3 [<K ij d t9ij - N l H g i\ 
— J dtNHgx + (boundary terms), (16) 
where 7r y is momentum conjugate to gij given by 

^ = T?^= M PlV9P ij , V' ffVW (17) 

o(Otgij) 

The Hamiltonian corresponding to the time t is the sum 
of constraints and boundary terms as 

H g [dt] = NH g± + J dx 3 N' l Hgi + (boundary terms). 

(18) 

Finally, by variation with respect to g^ (t, x), we obtain 
dynamical equation 



£gij "T" £rnij 0? 



(19) 



^9 ij — QikQjl 



2 SIg 

N V9 S 9ki ' 
2 SI m 



The explicit expression for E g ij is given by 
1 



(20) 



So 



Mpi 



1 



N 



(fit - N"D k ) Pi 



-(pikDjN" + p jk DiN") - K Pij + 2K« Pkj 



+ ^gijK kl p M + ^Agtj - Gi 



+ £z>lij, (21) 



where £ z >lij is the contribution from L z> \ and Gij is 
Einstein tensor of g^ . 

The invariance of I a under the infinitesimal transfor- 
mation §5§ leads to the following conservation equations, 
where a represents g or m. 



= Nd t H aJ _ 



dx A 



= -(dt-NW^Hai + KHci 



,(22) 



(23) 



B. Properties and problems of Hofava gravity 

Renormalization group (RG) flow of X 
The IR limit of the theory is characterized by the 
parameter A. When A = 1, the gravitational action 
in the IR limit is identical to the ADM form of the 
Einstein-Hilbert action except that the lapse function 
is independent of spatial coordinates. Note, how- 
ever, that setting A = 1 at all scales would lead to 
a problem since deviation from A = 1 in the UV is 
essential for avoidance of codimension-one caustics 
While the (RG) flow of Hofava-Lifshitz gravity has not 
yet been analyzed, we therefore suppose that A = 1 
is an IR fixed point of the RG flow so that A is suffi- 
ciently close to 1 in the IR but deviates from 1 in the UV. 

Instabilities 

In order to avoid ghost instability of the scalar graviton, 
the coefficient of the time kinetic term (A — 1)/(3A — 1) 
must be positive [l| . This makes the sound speed squared 
c 2 s = —(A — 1)/(3A — 1) negative [l3[ and leads to an IR 
instability in the linear level. As shown in Appendix [Cj 
this type of instability does not show up if 



c s \ < max HL, (ML) -1 , |$ 



1/2 



(24) 



where M is the energy scale suppressing higher-derivative 
terms, L (> 0.01mm) is the length scale of interest, H 
is the Hubble expansion rate at the time of interest and 
$ is the Newtonian potential at the position of interest. 
This is a condition on the way how A depends on L, 
H and $ and, thus, validity of this condition can in 
principle be checked by analyzing the RG flow. 

Strong coupling 
The strong coupling between the scalar graviton and 
matter pointed out in Q in the limit A — > 1 is absent if 
the projectability condition is imposed On the other 
hand, even with the projectability condition, there still 
remains the strong self-coupling of the scalar graviton 
in the limit A — > 1, as is clear from Hofava's original 
paper [l[. However, this is not a problem if the scalar 
graviton decouples from the rest of the world at the 
nonlinear level. In massive gravity theories flol | this kind 
of decoupling phenomenon due to nonlinear dynamics is 
known as Vainstein effect 111]]. A potential problem is 



that the strong self-coupling makes quantum corrections 
to the classical action very large. In non-renormalizable 
theories like massive gravity, this can be a fatal flaw 
since we really need to deal with an unknown quantum 
action including quantum gravity effects while the 
Vainstein effect has been shown only for a classical 
action with finite number of terms. On the other hand, 
in Hofava-Lifshitz gravity, if the theory is renormalizable 
then we can safely use the renormalizable action with 
finite number of terms to investigate whether the scalar 
graviton really decouples or not. Detailed investigation 
of decoupling of the scalar graviton is one of the most 
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important issues, but the present paper aims to shed 
light on a different aspect of the theory. 

Black hole 

With A = 1 and vanishing cosmological constant, 
Schwarzschild spacetime is locally an exact vacuum solu- 
tion of Hof ava-Lifshitz gravity (See appendix [A} . In the 
Painleve-Gullstrand coordinate system the Schwarzschild 
metric is 



ds 2 



-dt 2 



dr± 



2M 
r 



dt) +r 2 dfll, (25) 



and satisfies the projectability condition. Here, dQ,?, is 
the metric of the unit sphere. Since the spatial metric 



dr 2 



'dfl 2 , is flat, higher spatial curvature terms do not 



change the solution at all. 

As we note in Section HI in this theory a black hole is 
not "black" in the UV and we can see the central singu- 
larity of a Schwarzschild black hole if it exists. Since at 
the singular point the basic equations (Q~4|) and (fl~9f are 
not satisfied, a Schwarzschild metric must be modified 
near the center. For the Schwarzschild metric (|25p . the 
extrinsic curvature of constant-time hypersurfaces be- 
comes large near the center and the system enters the 
UV regime. Therefore, A should deviate from 1 and the 
Schwarzschild spacetime is no more a valid description 
of the geometry near the center. For this reason, the 
apparent singularity at the center is not physical. 



III. NO SPHERICALLY SYMMETRIC AND 
STATIC SOLUTION 

As the last topic of Section [ill we touched upon the 
problem of the central singularity in a black hole. In 
order to investigate what happens inside the black hole, 
we need to study gravitational collapse and formation of 
a black hole. For this purpose we need to formulate a 
regular initial condition. As a first step, in the present 
paper we consider a static star. Surprisingly enough, 
we find a no-go result: a spherically-symmetric, globally 
static solution can not be constructed in Hofava-Lifshitz 
gravity under the assumption that the energy density 
is a piecewise-continuous, non-negative function of the 
pressure and that the pressure at the center is positive. 
This section provides a proof of this statement. 



A. Painleve-Gullstrand coordinate 

We consider spherical-symmetric and static configura- 
tions. Since the lapse function does not depend on spatial 
coordinates, we can set it to unity by space-independent 
time reparametrization: 



For a spherically symmetric, static configuration, we can 
express the shift vector and the spatial metric as 

N i di = P(x)d x , g^dx'dx 1 = dx 2 + r 2 (x)dnl, (27) 

where dfl?, is the metric of the unit sphere and x is the 
proper distance from the center. Non-vanishing compo- 
nents of the Ricci tensor and Ricci scalar for the three- 
dimensional geometry are 

R% = R$=l[l-rr"-(r'f], (28) 



and the extrinsic curvature and its trace are 



(29) 



(r 2 BY 

Kijdx l dx 3 = -fldx 2 - f3rr'dn 2 2 , K = (30) 

where a prime denotes derivative w.r.t. x. The corre- 
sponding ADM metric is 

gffldx^dx" = -dt 2 + [dx + /3(x)dt} 2 + r 2 (x)dnj. (31) 

This is an analogue of the Painleve-Gullstrand coordinate 
system. Note that 



(32) 



is a timelike Killing vector. Global staticity requires £^ 
to be globally timelike, i.e. 1 — f3 2 > everywhere. The 
unit vector n M normal to the constant time hypersurface 
is given by 



n^dx^ = -dt, n^df, = d t - /3d x 



(33) 



It is instructive to see that the Newton potential 4>(r) 
and the Misner-Sharp energy m(r) are written as 



(34) 



= g^{dr),{dr) v = {l-p 2 ){r'f-{^) 



B. Matter sector 

As for (real) matter, for simplicity we consider the 
perfect-fluid form which is at rest w.r.t. the Killing vec- 
tor 



T MV = p{x)u^u v + P(x) 

11 _ s 



w = 



(36) 
(37) 



N = l. 



(26) 



Its components relevant for the ADM decomposition are 

T^n v = (l-j3 2 )-\p + P)-P, (38) 
T^dx 1 = (1- p 2 )- l (p + P)/3dx, (39) 
T l3 dx l dx j = [P 2 p + P]dx 2 + Pr 2 dVL\. (40) 
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C. Inconsistency near the center 

In this subsection, we show that a spherically- 
symmetric compact object cannot be globally-static un- 
der a set of reasonable assumptions in the matter sec- 
tor. We assume that the energy density p is a piecewise- 
continuous 1 non-negative function of the pressure P and 
that the pressure at the center P c is positive. The three- 
dimensional spatial geometry and the extrinsic curvature 
must be regular at the center because the constant-time 
surfaces are physically embeded in Hofava-Lifshitz grav- 
ity. (See Appendix IDl on this point.) We prove that the 
regularity condition at the center is incompatible with 
the momentum conservation law of the matter under the 
above assumptions. Since the momentum conservation 
equation does not include higher spatial curvature terms 
in the gravity action, our proof is totally insensitive to 
the structure of higher spatial curvature terms and holds 
for any values of the dynamical critical exponent z. 

In the following argument, we make the proposition 
that there exists a spherically-symmetric, globally-static, 
regular solution, and show contradiction. As commented 
after eq. (|32|) . the global-staticity implies that (1 — (P) is 
positive everywhere. 

The momentum conservation equation, (I23p with a = 
m, becomes 



P'(l-/3 2 ) + (p + P)(l-/3 2 )' = 0. 



(41) 



The regularity of the extrinsic curvature (I30[) implies that 
j3' is finite. This and (|4"Tj) imply that P' is also finite. As 
a corollary, /3 and P are continuous functions of x. Since 
p is assumed to be a piecewise continuous function of P, 
this means that p + P is a piecewise continuous function 
of x. 

Let x c be the value of x at the center. Since we have 
assumed that p is non-negative everywhere and that P c > 
0, the continuity of P(x) implies that p + P is positive 
in a neighborhood of the center. Now let us define xq as 
the minimal value for which at least one of (p + P)\ x=xo , 
lim x ^ Xo ^ (p + P) and lim x _>. Xo+ o (p + P) is non-positive. 

Dividing eq. (|4ip by (p + P)(l — /3 2 ) and integrating it 
over the interval x c < x < xq, we obtain 



ln(l - fil) - ln(l - Pi) 



x -0 



P' 



p + P 



dx, (42) 



where j3 c = /3(x — x c ) and /3q = (3{x = xq). The regu- 
larity of the Ricci scalar (|2U)) and the extrinsic curvature 
(|30|) at the center implies that r' c = 1 and f3 c = 0, where 
r' c is the value of r' at the center. Therefore, the left 
hand side of eq. (|42|) is non-positive. 



We assume piecewise-continuity instead of continuity, in order 
to allow dP/dp to vanish in a finite interval of p. 



Since P is a differentiable function of x, the right hand 
side of eq. (|4"2"j) can be transformed as 



Xn-0 



P' 



P + P 



-dx 



dP 



P c P(P) + P' 



(43) 



where P = P(x — Xq). The definition of xq implies 
that at least one of (p + P)\ x =x i ^ m x~^xa-o{p + P) an d 
lim x ^ Xo+ o{p+ P) is non-positive. Since we have assumed 
that p is non- negative everywhere, Pq = lirn x _>a; _o P = 
linix-^xo+O P is non-positive. Thus, we have 



Po < < P c 



(44) 



This implies positivity of the right hand side of (|4"5|) since 
from the definition of Xo the integrand is positive in the 
domain of integration. This leads to an contradiction 
with the previous statement that the left hand side of 
(|4"2"|) should be non-positive. 



IV. SUMMARY AND DISCUSSION 

In Hofava-Lifshitz gravity, regularity of a solution re- 
quires smoothness of not only the spacetime geometry 
but also the foliation. As a result, the regularity con- 
dition at the center of a star is more restrictive than in 
general relativity. Under the assumptions that the energy 
density is piecewise-continuous non-negative function of 
the pressure and that the pressure at the center is pos- 
itive, we have proved that the momentum conservation 
law is incompatible with the regularity at the center for 
any spherically-symmetric, globally-static configurations. 
The proof is totally insensitive to the structure of higher 
spatial curvature terms and, thus, holds for any values 
of the dynamical critical exponent z. Therefore, under 
the assumption we made on the matter sector, we con- 
clude that a spherically-symmetric star should include a 
time-dependent region, presumably near the center. 

The assumptions we made are physically natural. For 
example, a polytropic fluid satisfies them. Note that our 
proof does not assume asymptotic flatness and that a 
cosmological constant A can be included in the gravity 
action ([7]). Shifting p, —P and —Mp t A with the same 
amount does not change the physical system but may 
validate/invalidate some of the assumptions of the proof. 
In order to construct a static star solution, we need to 
violate at least one of the assumptions for all possible 
choices of such a shift. One possibility is to introduce an 
exotic matter such as a quintessence field. Introduction 
of an exotic matter is, however, not necessarily sufficient 
for the existence of a static star solution. 

One must not consider our result, i.e. nonexistence of 
static star, as a serious problem of Hofava-Lifshitz the- 
ory. It is known that this theory does not recover general 
relativity at low energy but can instead mimic general 



6 



relativity plus cold dark matter ji[ 2 The existence of 
built-in "cold dark matter" is an inevitable prediction of 
the theory and might solve the mystery of dark matter 
in the universe. Our result in the present paper is to- 
tally consistent with this picture: as in the standard cold 
dark matter scenario, the "cold dark matter" accretes 
toward a star and thus inevitably makes the stellar cen- 
ter dynamical. This is the physical reason why there is 
no static star in Hof ava-Lifshitz theory and, thus, our re- 
sult strongly supports the "dark matter as an integration 
constant" scenario Q. 

In subsection III BL we have commented on the Vain- 
stein effect for the scalar graviton, and also derived the 
condition (|24l) under which linear instability of the scalar 
graviton does not show up. This condition should be con- 
sidered as a phenomenological constraint on properties of 
the renormalization group (RG) flow. Detailed analysis 
of the RG flow and the fate of scalar graviton will be one 
of important future subjects. 



Acknowledgments 

The authors would like to thank K. Takahashi for 
useful discussions. The work of S.M. was supported 
in part by Grant-in-Aid for Young Scientists (B) 
No. 17740134, Grant-in-Aid for Creative Scientific Re- 
search No. 19GS0219, Grant-in-Aid for Scientific Re- 
search on Innovative Areas No. 21111006, Grant-in-Aid 
for Scientific Research (C) No. 21540278, and the Mit- 
subishi Foundation. This work was supported by World 
Premier International Research Center Initiative (WPI 
Initiative), MEXT, Japan. 



Appendix A: Asymptotically-flat vacuum solution 
with A = 1 and without cosmological constant 

For T^„ — 0, A — and A = 1, the momentum con- 
straint (Ti"4l and the dynamical equation (fT9"|) become 

= p—, (Al) 
r 

0=i [2/3rr"-(l-/3 2 )(r') 2 

+2/3/3Vr' + l] + i(r') 2 ^ 1 , (A2) 

where ££>i = 9 ik 9 il S z> iki- 

On the other hand, we shall not impose the Hamilto- 
nian constraint equation ([9]) to our solutions. The reason 
is as follows. In Hofava-Lifshitz gravity the Hamiltonian 



2 The constraint algebra is smaller than in general relativity since 
the time slicing is synchronized with the rest frame of cold dark 
matter in the theory level. 



constraint is not a local equation satisfied at each spatial 
point but an equation integrated over a whole space. On 
the other hand, our assumption that — is justified 
only if we consider physics at length scales sufficiently 
shorter than the cosmological horizon. Otherwise, the 
stress-energy tensor of cosmological fluids should be in- 
cluded. Also, our assumption of staticity completely ne- 
glects cosmological expansion of our universe. Therefore, 
static vacuum solutions would be valid only in a (large 
but) finite region. For this reason, it is not appropri- 
ate to substitute a static vacuum solution to the global 
Hamiltonian constraint, which is an integral over a whole 
space including the region outside the cosmological hori- 
zon. The contribution from the region within the regime 
of validity of a static vacuum solution to the integral can 
easily be canceled by the contribution from the region 
outside the regime of validity of the solution. Therefore, 
we must not impose the Hamiltonian constraint to our 
static vacuum solution. 

There are two branches of solutions to the momentum 
constraint equation (|A1|) : j3 = and r' = r\ (constant). 
In the following we shall consider each branch separately. 
Without loss of generality, we assume that n > Q. (If 
n < then we can make redefinition x — > — x so that 
ri -n.) 

Let us consider the first branch, /3 = Q. In this case 
the dynamical equation (| A2|) is reduced to 

0=i[l-(r') 2 ] + 5(r') 2 ^i. (A3) 

For asymptotically-flat solutions, the first term behaves 
as ~ 1 /r 2 unless r' = ± 1 , while £ z decays at least as fast 
as 1/r 4 in the asymptotic region. Therefore, asymptot- 
ically flat solutions with (3 — has r' = ±1 and can be 
included in the second branch of solutions to (|Al|) . i.e. 
r' = t*i (constant). 

In the second branch, r' = r\ (constant), the Ricci 
tensor and the Ricci scalar become 

1 — t- 2 2 — 2r 2 

R% = 0, R§ = ^fi, R=^f 1 - (A4) 

Evidently, the value n = 1 is special. In this case, all 
components of the Ricci tensor vanish and, thus, £^>i = 
0. Thus, for n = 1, (|A"2j) is reduced to (r/3 2 )' = and 
leads to the Schwarzschild solution (|2"5|) . 

For r' — r\ (constant) with r± =/= 1, we can see from 
eq. (|A4[) that those terms in £f>i including 2n spatial 
derivatives are proportional to 1/r 2 ™, where n = 2, • • • ,z. 
Thus, the dynamical equation (|A2[) is written as 

r^r)'=rl-l-±^, (A5) 

where a n {r\) (n = 2, • • • , z) are constants depending on 
the constant r\ and the structure of the higher spatial 
curvature terms I z> \. Since all components of the Ricci 
tensor vanish for r% = 1, we have 

a„(l)=0, (n = 2,---,z). (A6) 
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Integrating the equation (|A5I) . we obtain 



2// 



2n 

n=2 



(A7) 



where /i is an integration constant. 

Note that /3 must be real and thus /3 2 > 0. Since (jA"7|) 



implies that /3 



1)/ r\ (r —> oo), it follows that 



ri > 1. (Note that we have assumed that r\ > without 
loss of generality.) 

In summary, the asymptotically-flat vacuum solution 
is characterized by two integration constants r\ (> 1) and 
p. The solution is given by r' = r*i and (|A7|) . When r\ = 
1, the solution is reduced to the Schwarzschild spacetime 
([25jl with Af = p since a„(l) vanishes as shown in (|A6|) , 



Appendix B: Yet another no-go result 

In this appendix we show yet another no-go result, us- 
ing the momentum constraint equation. As in scction lHll 
we make the proposition that there exists a spherically- 
symmetric, globally-static, regular solution under a set 
of assumptions, and show contradiction. We set A = 1 
and restrict our consideration to asymptotically flat so- 
lutions (thus with vanishing cosmological constant). We 
also assume that the shift j3 is analytic. As for the matter 
sector, we only assume the null energy condition. 

The momentum constraint equation is 



r" + 



p + P 



2M 2 pi 1 - /3 2 



= 



(Bl) 



and is insensitive to the structure of higher spatial cur- 
vature terms. If /3 = in a finite interval of x then /3 = 
everywhere since /3 is assumed to be analytic. This cor- 
responds to a trivial flat solution. We therefore suppose 
that 



P 



P 



2M pl 1 



/? 2 







(B2) 



is satisfied everywhere. If the null energy condition (p + 
P > 0) is satisfied then by integrating (|B2[) from the 
center towards outside we obtain 



P- 



P 



2M pl 1 



dx < 0, (B3) 



where X — X q IS the center and x ou t > x c . The equality 
holds if and only if p + P = everywhere in the region 
x c < x < x ou t- Note that, as already stated just after 
(|52")l . the global staticity requires 1 — (3 2 > 0. Since we are 
interested in a stellar solution, we suppose that p + P =/= 
somewhere. Therefore, we obtain 



< r' 



for sufficiently large x ou t- (B4) 



Now we can show that the regularity of the solution at 
the center is incompatible with the asymptotic flatness. 



The regularity at the center requires that r'\ x=Xc = 1. 
On the other hand, in Appendix [X] we see that the 
asymptotically-flat vacuum solution always has r' > 
1. Therefore, the stellar solution can be connected to 
the vacuum solution in the asymptotic region only if 
r'\ x =x out > 1) where we have chosen a sufficiently large 
x ou t. This is in conflict with (|B4[) . 



Appendix C: Stealth linear instability of scalar 
graviton 

As shown in j~0 • the lack of local Hamiltonian con- 
straint leads to "dark matter as an integration constant" , 
a non-dynamical component which behaves like pressure- 
less dust. As in the standard CDM scenario, the dust- 
like component exhibits Jeans instability and forms large- 
scale structures in the universe. The timescale of Jeans 
instability is 



t, 



(Gnp) 



-1/2 



(CI) 



where Gjv is Newton constant and p is the energy density 
at the position of interest. Note that this instability is 
necessary for structure formation if we consider the dust- 
like component as an alternative to CDM. 

As mentioned in the introduction, the scalar graviton 
exhibits linear instability due to the negative sound speed 
squared, c 2 = —(A — 1)/(3A— 1) < 0. The corresponding 
time scale is 



L 



(C2) 



where L is the length scale of interest. Thus, as far as 

t L > tj, (C3) 

the linear instability does not show up. 
Since the dispersion relation is of the form 

oj 2 = k 2 x [c 2 + 0(k 2 /M 2 )] , (C4) 

the linear instability is stabilized by higher derivative 
terms if 



< 



1 



ML' 



(C5) 



where M is the energy scale suppressing higher derivative 
terms and we have assumed that the coefficient of k 2 /M 2 
in the square bracket in (|C4|1 is positive and of order 
unity. Also, the linear instability is tamed by Hubble 
friction if 



(C6) 



where H is the Hubble expansion rate at the time of 
interest. 

In summary, if one or more of the three conditions 
(|C3[) , (|C5|) and (|C6[) is satisfied then the linear instability 
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of the scalar graviton does not show up. By introducing 
Newton potential $ ~ —GpjpL 2 , these conditions are 
summarized as 



< max 



HL, {ML)- 1 , |$ 



1/2 



(C7) 



As argued in [4], nonlinear extension of the linear in- 
stability should be formation of would-be caustics and 
bounce. For length scales shorter than ~ 0.01mm, we 
do not experimentally know how gravity behaves and, 
thus, the existence of formation of would-be caustics and 
bounce does not contradict with any experiments. In 
some sense, this is similar to the so called spacetime foam. 
Therefore, in (|C7|) . we do not have to consider the length 
scale L shorter than ~ 0.01mm. 



Appendix D: Non-triviality of regular foliation 



Throughout this paper (eg. in subsection IIII CI) , we 
demand the regularity of time foliation as a physical con- 
dition. In this appendix we shall see how significant this 
condition is. For this purpose, we consider a simple stel- 
lar solution in general relativity in Painleve-Gullstrand- 
like coordinate system and show that the time foliation is 
actually singular at the center. This makes it clear that 
the regularity of time foliation is not a trivial but signifi- 
cant condition. Based on this observation, we shall fore- 
see what happens near a stellar center in Hofava-Lifshitz 
gravity. 

One of the simplest (and idealized) stellar solutions 
in general relativity is that with uniform energy den- 
sity. The explicit expression can be found in textbooks, 
e.g. [l2j]- By going to a Painleve-Gullstrand-like coordi- 
nate system, the solution is written as 



ds 2 

P 
P 
Po 



-dt 2 + e~ 2lf '{dr + Pdtf + r 2 dnl, (Dl) 
po = constant, (D2) 
(1 - 2M/R) 1 / 2 - (1 - 2Mr 2 /R 3 ) 1 / 2 



(1 - 2Mr 2 /i? 3 ) 1/2 - 3(1 - 2M/R) 1 / 2 



(D3) 



where 



&i°s(-K 1 ->/5^)) ( r<jB W 

(r>R) 



^-r 2 Pa + - 1 (r < R) 



(r > R) 



47T 

M = —R 3 p , 



(D5) 
(D6) 



and r — R is the surface of the star. For this solution, 
the Ricci scalar of the three-dimensional spatial metric is 



R= %{l-e-^ +2e~^A, (D7) 



and diverges at the center. It is also easy to see that 
the extrinsic curvature of constant-time hypersurfaces 
also diverges at the center. Physical reason for these 
divergences is easy to understand. A metric with a con- 
stant lapse function is characterized by a congruence of 
geodesies orthogonal to constant-time hypersurfaces. In 
general relativity, a contracting congruence of geodesies 
forms caustics. This is the physical reason why the three- 
dimensional spatial geometry and the extrinsic curvature 
are singular at the center. 



In general relativity, divergence of the extrinsic cur- 
vature and the three-dimensional spatial curvature can 
be just a coordinate singularity. Indeed, one can easily 
see that the four-dimensional geometry is regular at the 
center of the above solution. Therefore, divergence at 
the center is just a coordinate singularity and is not a 
problem for the above solution in general relativity. 



On the other hand, in Hofava-Lifshitz gravity, diver- 
gence of the three-dimensional spatial curvature or/and 
the extrinsic curvature is a physical singularity. Thus, 
the above solution with constant energy density is not 
physically viable. Indeed, since the three-dimensional 
spatial curvature diverges, higher spatial curvature terms 
become important near the center and should change the 
behavior of the solution. It is also expected that devia- 
tion of A from 1 should also be important near the center. 
Those effects should generate "dark matter as integration 
constant" , which should develop would-be caustics and 
bounce [4j. What we expect is occurring near the center 
is, thus, a sequence of microscopic would-be caustics and 
bounces. 
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